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We examine the interplay of nonlinearity of a dynamical system and thermal fluctuation of its 
environment in the "physical limit" of small damping and slow diffusion in a semiclassical context 
and show that the trajectories of c-number variables exhibit dynamical chaos due to the thermal 
fluctuations of the bath. 
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PACS number(s): 05.45.Mt 

The interplay of nonlinearity of a dynamical system 
and thermal fluctuations of its environment has been 
one of the major areas of investigation in the recent 
past. These studies have enriched our understanding 
of nonequilibrium processes in several contexts, such as, 
symmetry between the growth and decay of classical fluc- 
tuations in equilibrium [EJ, interesting topological fea- 
tures of patterns of paths of large fluctuations in nonlin- 
ear systems [El , existence of generalized nonequilibrium 
potential influence of nonlinearity on dissipation in 
multiphoton processes H and higher order diffusion in 
a nonlinear system etc. While the development in 
these areas is largely confined to classical domain we ex- 
amine a related issue in the semiclassical context. Since 
quantization is likely to add a new dimension to the in- 
terplay of nonlinearity and stochasticity in a weakly dis- 
sipative system, it is worthwhile to consider the physical 
limit of small damping and slow diffusion due to thermal 
fluctuations of the environment and look for the ther- 
mal fluctuations-induced features of nonlinearity in the 
dynamics. In this communication we specifically explore 
some interesting aspects of dynamical chaos in a driven 
bistable system whose origin lies at the fluctuations of 
the environmental degrees of freedom. 

To describe the dissipative quantum dynamics of a 
system we first consider the traditional system-reservoir 
model developed over the last few decades ||,|5|||- The 
Hamiltonian of the bare system is coupled to an environ- 
ment modeled by a reservoir of harmonic oscillator modes 
characterized by a frequency set {fij}. The quantum dy- 
namics is generated by the overall Hamiltonian operator 
H for the system, environment and their coupling as fol- 
lows; 
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defines the usual kinetic and potential energy terms cor- 
responding to the system, x and p being the position and 



momentum operators, respectively. The second and the 
third terms in (1) specify the reservoir modes and their 
linear coupling to the system. g(£l) denotes the system- 
reservoir coupling constant. 

Systematic elimination of the reservoir modes in the 
usual way, using Born and Markov approximations leads 
one to the following standard reduced density matrix 
equation of the system only Q . 
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Here the system operator co-ordinate x is related to the 
creation and the annihilation operators a\ a respectively 
as x — (-^i=)(a + a^). a; is the linearised frequency of 
the system. Also the spectral density function of the 
reservoir is replaced by a continuous density T>(lo). 7 > 



is the limit of 27t\g(cu 



as u> —> + and is assumed 



to be finite. 7 is the relaxation or dissipation rate and 
nj is the diffusion coefficient D. h(— [exp (jpf) — 
is the average thermal photon number of the reservoir. 
The terms analogous to Stark and Lamb shifts have been 
neglected. 

The first term in Eq.(3) corresponding to the dynam- 
ical motion of the system refers to Liouville flow. The 
terms containing 7 arise due to the interaction of the sys- 
tem with the surroundings. While ^ terms denote the 
diffusion of fluctuations of the reservoir modes into the 
system mode, ^ terms refer to the loss of energy from the 
system into reservoir. In the limit T — > 0, i.e., tt, — > the 
system is influenced by pure quantum noise or vacuum 
fluctuations [ Note that 1 in (n+1) in Eq.(3) corresponds 
to the vacuum] . 

Our next task is to go over from a full quantum op- 
erator problem described by the Eq.(3) to an equivalent 
c-number problem . To this end we consider the quasi- 
classical distribution function W (x,p, t) of Wigner. x,p 
are now c-number variables. Rewriting Eq.(3) in a quasi- 
classical language we obtain 
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77 in Eq.7 is a parameter used to identify the environment- 
induced effect on the dynamics described by Eq.(4) (kept 
for bookkeeping in the calculation and put rj = 1 at the 
end). 

In the semiclassical limit Tilj « kT, we have n + 1 ^ h 
and D « 'ykT so that Eq.(7) reduces to 
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The overall dynamics described above is a superposi- 
tion of two contributions, i. e. , the Liouville-Wigner 
dynamics and the system-reservoir dissipative dynam- 
ics. That the two contributions act independently is 
an assumption. The master equation (3) [or its Wigner 
function version (4)] is the most popular one in quan- 
tum optics. It has been extensively used 0, for the 
strongly nonlinear processes like three-wave, four-wave 
mixing and strong coherent light-matter interaction phe- 
nomena. The equation has also been applied in the con- 
text of chaos, e. g., in the dissipative standard map 
dissipative logistic map || , semiclassical theory of quan- 
tum noise in open chaotic systems and in the 
studies of decoherence in relation to chaos for analysis 
of quantum-classical correspondence Jl2|,|l3| . 

In the semiclassical {% — > 0) limit the dissipative quan- 
tum dynamics can be conveniently described by "WKB- 
like" ansatz (we refer to " WKB-like" since we are consid- 
ering more that one dimension. Traditional WKB refer 
to one dimension only) Jj],[l4| of Eq.(5) for Wigner func- 
tion of the form 



W(x,p,t)=Z(t)exp(--) 
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Here Z(t) is a prefactor and s(x,p, t) is a classical ac- 
tion which is a function of c-number variables x and p , 
satisfying the following Hamilton-Jacobi equation 
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The derivation of Eq.(7) is based on the consideration 
of the "physical limit" of weak dissipation and slow fluc- 

(note that 



tuations in the sense §r 



j- where D\ = 



D\ and h have same dimension). The above equation 
can be solved by integrating the Hamiltonian equations 
of motion, 
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which are derived from the following Hamiltonian H e f f 
H ef f =pP- V'X - ^7](xP +pX) + f](P 2 + oj 2 X 2 ) 



E(_1\3n+1 Q2n+lT/ 
\ l ) £ _X 2n+1 
2 2 "(2n+ 1)! dx 2 '^ 1 

n>l 



(9) 



Here we have put §^ = P and |p = X. The introduction 
of additional degree-of- freedom by incorporating the aux- 
iliary momentum (P) the and co-ordinate (X) makes the 
system effectively a two-degree-of-freedom system. The 
origin of these two variables is the thermal fluctuations 
of the reservoir. It is easy to identify the environment- 
related terms containing rj in Eqs.(7-9). The auxiliary 
Hamiltonian (9) is therefore not to be confused with the 
microscopic Hamiltonian (1) which describes a system in 
contact with a reservoir with infinite degrees of freedom. 
Although the phase space trajectories concern fluctua- 
tions of c-number variables in the formal sense, because 
of the equations of motion (8) described by a Hamilto- 
nian (9), the motion is strictly deterministic. The ex- 
periments on the corresponding classical version of the 
problem by Luchinsky and McClintock jj| have demon- 
strated that a trajectory of fluctuation is indeed a part 
of physical reality. We emphasize, however, here a num- 
ber of distinguishing features in this context. While the 
studies by Luchinsky and McClintock jj] and Graham 
and Tel Q concern overdamped limit, we consider here 
a weakly dissipative system. Furthermore because of the 
quantum correction, the phase space trajectories of fluc- 
tuations are significantly modified by semiclassical fea- 
tures. The introduction of these quantum features at a 
semiclassical level through a c-number Hamiltonian de- 
scription of a dissipative evolution in the physical limit 
of weak damping and slow diffusion due to thermal noise 
is the essential content of the present work. Since the 
Eq.(9) describes deterministic evolution under nonlinear 
potential, the pattern of trajectories of fluctuations may 
display chaotic behaviour. In what follows we investigate 
this dynamical aspect of the dissipative system. 

The equations derived in the weak thermal noise limit 
for the weakly dissipative semiclassical systems are fairly 
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general. For illustration we now consider a simple model 
system Hamiltonian Ho (see Eq.2) 

2 

Hq = — + ax* — bx 2 + gxcosQot (10) 

which describes a bistable potential driven by a time- 
periodic held, a and b are the constants of the potential 
V(x). The fourth term in (10) includes the effect of cou- 
pling of the system as well as the strength of the held of 
frequency f2o- For the Hamiltonian (10) the Eqs.(8) read 
as 

1 

x = p — —rjx 2?]P 

X = P - -riX 
2 ' 

p = 4ax 3 -2bx + g cos fit + r](^p - 2uj 2 X) - 3axX 2 
P = (12ax 2 - 2b)X + ^tjP - aX 3 (11) 

which are derivable from the auxiliary Hamiltonian 

H ef f = Pp- (4ax 3 - 2bx + gcosujt)X - ^r](xP + pX) 
+r](P 2 + LU 2 X 2 ) + axX 3 . (12) 

The system Hamiltonian (10) has served as a standard 
paradigm for a number of theoretical and experimental 
investigations |l],|l0|,|l5| over many years. For the present 
purpose we choose the following parameter values; a = j , 
b = |, fl — 6.07 and g = 10. Since we are considering 
the physical limit of weak damping and small diffusion 
we take the initial conditions for the auxiliary variables 
X and P (which originate from the fluctuations due to 
environment) as P = and X — > (we have used X = 
1.5 x 10~ 6 ). This ensures a vanishing Hamiltonian H e ff 
for the entire numerical investigation that follows below. 

We hrst consider a specific trajectory with the initial 
condition p = and x — —2.512 for small values of 7 
(typical 0.01). Under this condition^ = 1) the system 
is vanishingly coupled to the surroundings and conse- 
quently the dynamical behaviour is effectively due to the 
weak dissipation only. We illustrate this situation in Fig. 
1 in terms of a Poincare map for the phase space which 
exhibits strong global chaos. On the other hand when 
the parameter 77 is switched off (77 = 0) the system dis- 
plays typical weak chaos (Fig. 2). The similar behaviour 
has been observed for other sets of initial conditions for 
x and p (We have not reproduced them here for the sake 
of brevity) . 

The effect of weak dissipation and slow diffusion due to 
thermal fluctuations from the surroundings can be seen 
in the case of other sets of initial conditions also. For the 
initial condition p = and x — —2.49 one observes for 
7 = 0.01 dissipative strong chaos (77 = 1). This is illus- 
trated in Fig. 3. It is interesting to note that for rj = the 
same trajectory gets localized in the left well as a regular 



one as shown in Fig. 4. The weakly chaotic and regu- 
lar trajectories in Figs. 2 and 4, respectively, are purely 
semiclassical in nature (in the absence of any coupling to 
the surroundings). The strong global chaotic behaviour 
as shown in Figs. 1 and 3 has therefore its origin in 
the thermal fluctuations of the reservoir. In other words 
the chaotic behaviour or its enhancement is exclusively 
due to thermal fluctuations from the surroundings, which 
becomes appreciable even in the physical limit of weak 
damping and slow diffusion of the thermal noise within 
the semiclassical description. We have checked this as- 
sertion for some other values of the initial conditions for 
the system oscillator. 

The reduction of the system-reservoir Hamiltonian de- 
scription [H in Eq.l] for a dissipative quantum system to 
an auxiliary Hamiltonian description [H e ff in Eq.9] effec- 
tively reduces the infinite-degree-of-freedom system to a 
two-degree-of-freedom system where the auxiliary degree 
of freedom characterized by X and P owe their origin in 
the fluctuations of the reservoir. Since X and P appear 
as the multiplicative factors in the auxiliary Hamiltonian 
H e f f , the weak thermal noise limit makes H e f t a vanish- 
ing Hamiltonian. The observed semiclassical chaos may 
therefore be regarded as a dynamical manifestation of the 
interplay of nonlinearity and thermal fluctuations. 

In this paper we have examined the weak thermal noise 
limit of a semiclassical dissipative nonlinear system. We 
have shown that the vanishing Hamiltonian method can 
be suitably extended to follow the phase space trajecto- 
ries of fluctuations of c-number variables which exhibit 
dynamical chaos. In view of the accessibility of the model 
to analogue electronic circuits fl we believe that the re- 
sults discussed bear further experimental relevance in the 
semiclassical context. 
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support and to J. Ray Chaudhuri for helpful discussions. 
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FIG. 1. Plot of x vs p on the Poincare surface of section 
(X = 0) for Eq.(ll) with initial condition x = -2.512, p = 0, 
X -> 0, P = 0, 7 = 0.01 r) = 1. (Units are arbitrary). 



FIG. 2. Same as in Fig.l but for n = 0. 



FIG. 3. Same as in Fig.l but for x = -2.49. 



FIG. 4. Same as in Fig.3 but for r\ = 0. 
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